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Extragalactic Cepheids are the basic rungs of the cosmic distance scale. They are excellent standard candles, although their luminosities 
and corresponding distance estimates can be affected by the particular properties of the host galaxy. Therefore, the accurate analysis of the 
Cepheid population in other galaxies, and notably in the Andromeda Galaxy (M31), is crucial to obtaining reliable distance determinations. 
We obtained accurate photometry (in B and V passbands) of 416 Cepheids in M31 over a five year campaign within a survey aimed at the 
detection of eclipsing binaries. The resulting Cepheid sample is the most complete in M 31 and has almost the same period distribution as the 
David Dunlap Observatory sample in the Milky Way. The large number of epochs (~ 250 per filter) has permitted the characterisation of the 
pulsation modes of 356 Cepheids, with 281 of them pulsating in the fundamental mode and 75 in the first overtone. The period-luminosity 
relationship of the fundamental mode Cepheids has been studied and a new approach has been used to estimate the effect of blending. We find 
that the blending contribution is as important as the metallicity correction when computing Cepheid distance determinations to M31 (~ 0.1 
mag). Since large amplitude Cepheids are less affected by blending, we have used those with an amplitude > 0.8 mag to derive a distance 
to M 31 of (m - M) = 24.32 ± 0.12 mag. 

Key words. Cepheids - stars: distances - galaxies: individual: M 3 1 - galaxies: distances and redshifts - methods:observational 



1. Introduction 

Cepheids are probably the most studied variable stars. Their 
large amplitudes and intrinsic luminosities make them eas- 
ily detectable in most photometric variability surveys. In ad- 
dition, their well-known p eriod-luminosity (P-L) relationship 
( Leavitt & Pickerindl912l) has made these variable stars one of 
the main cornerstone s in deriving extragalactic distances (see 
Freedman et al 2001, for an historical review). The importance 



1929h . After the observations of iBaade & Swopd dl965l and 
references therein), few efforts have been dedicated to further 
analyse the Cepheid population in M 3 1 . 

This trend is changing in recent years with the emer- 
gence of new observational capabilities. Several variabilit y sur- 
veys have started to study the stellar content in M 3 1 (IMacri 
2004[ and referen ces th erein) and other Local Group galax 



of Cepheids for distance determination stands in contrast with 
the relative lack of additional information on the specific char- 
acteristics of extragalactic Cepheids and the possible correc- 
tions because of their particular properties (i.e., metallicity). A 
clear example is the Andromeda galaxy (M31), where the first 
identification of Cepheids was already performed by Hubble 



ies (U dalski et al. 1999; Mocheiska et al. 2001- JlJdalski et al. 
200 ll: IPietrzvriski et alJ |2004 IMacri et alj l200fj|) . obtaining 



* Based on observations made with the Isaac Newton Telescope op- 
erated on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astroffsica de Canarias. 



large samples of Cepheids with accurate photometry. The de- 
tailed study of the observed Cepheids has emphasized the im- 
portance of an issue that was usually overlooked in most pho- 
tometric studies: the ef fect of bl ending. It has been proposed 
dMocheiska et al.l200fi hereafter lMMSS OO) that the magnitude 
of Cepheids may be affected by the light of unresolved com- 
panion stars (i.e., blends). The effect of blending is somewhat 
different from crowding or confusion noise, since companion 
stars appear to be in the same point-like source. Therefore, 
even when achieving a perfect point-spread function model- 



2 



F. Vilardell et aL: A comprehensive study of Cepheid variables in the Andromeda galaxy 



ing, blending could still be present. The effect can be the same 
as in spectroscopic binaries, where the individual components 
cannot usually be resolved from ground-based images. 

When studying extragalactic Cepheids the spatial resolu- 
tion decreases linearly with the distance of the host galaxy and, 
as a consequence, the number of possible blends increases. 
Small blending values are expected for the Large Magellanic 
Cloud, where individual Cepheids can be resolved from neigh- 
boring stars, even from ground-based observations. The situ- 
ation changes in M31 and M33, where mean blending val- 
ues of 0.18 mag and 0.16 m ag, respec tively, hav e been ob- 
tained dMocheiska et al.ll2004 hereafter ImMSS04 . These re- 
sults would imply, when extrapolated to more distant galax- 
ies, a downward revision of the Hubble constan t between 
5% to 10% (as explained by iGibson et all |2000|) . By con- 



trast, iBresolin et al.l d2005l) found an upper limit on blend- 
ing in NGC300 ( at ~ 2 Mpc) of 0.04 mag. In addition, 



(at 
20)! 



Gibs on et al.l (12000) showed that the systematic effect of blend- 
ing on farther galaxies (between 4 and 25 Mpc) is almost neg- 
ligible. Therefore, results obtained so far seem to indicate that 
blending has an important contribution when observing Local 
Gro up galaxies and dim inishes when observing distant galax- 



ies. IGibson et al.l yOOO) explained such behavior as a conse- 



quence of the background levels in M31 and M33 (with a 
large number of stars detected around Cepheid variables) not 
being representative of the more distant galaxies (where the 
background levels are the result of several unresolved sources). 
Because of the importance of the subject, a comprehensive 
study of the effect of blending on Cepheid distance determi- 
nations is highly desirable. 

One of the most recent variability surveys in M 3 1 was car- 
ried out by us an d was centered on th e Norm-Eastern quad- 
rant of the galaxy ( Vila rdell et al J2OO6I hereafter PaperJ)- Our 
main goal was to detect eclipsing binaries, although a sample 
of 416 Cepheids also resulted from the reduction and analysis 
process. The large number of detected stars and the high qual- 
ity of the resulting light curves have allowed the current study 
of Cepheid properties in M 3 1 . 



2. Photometric sample 

The photometric survey in Paper I was carried out with the 



2.5 m Isaac Newton Telescope (INT) at La Palma (Spain) 
over the course of five observing runs (between 1999 and 
2003). An implementatio n of the difference image analysis 
technique (Wozniak 2000) was used in the reduction process 
to automatically detect 3 964 variable stars with ~250 obser- 
vations _gej_band_(_B_jnd_J£^TJieanalysis of variance algo- 
rithm ( Schwarzenberg- Czernvll 19961) was used to identify pe- 
riodic variable stars. The folded light curves were visually in- 
spected and lead to the identification of 437 eclipsing bina- 
ries and 416 Cepheids. The photometric measurements were 
transform ed to the stand ard system through the observation 
of several iLandolj (119921) star fields, obtaining magnitude er- 
rors between 0.03 and 0.11 mag. Final ly, the phase-weighted 
mean magnitudes (ISaha & Hoesselll 19901) of the Cepheid sam- 
ple were determined for each passband (B and V). Figure Q] 
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Fig. 1. Color-magnitude diagram for 37 241 objects with photo- 
metric errors below 0. 1 mag in the reference catalog of Paper J 
(grey dots) and the 416 identified Cepheids (black circles). 
Objects above the dashed line are saturated. Error bars on the 
left indicate the mean error, in magnitude and color, at different 
magnitude values. 



shows the color-magnitude diagram of the observed stars with 
good photometric precision. Cepheid variables are highlighted. 

3. Period distribution 

The period distribution of Cepheids in Local Group galax- 
ies has been a major source of debate. Evolution models pre- 
dict a displacement on the peak of the period d istribution as 



a fun ction of the metallicity of the host galaxy (Bec ker et al 



19771) . It has also been observed that the Milky Way (MW) pe- 



riod distribution displays a dip at around 10 days, while such 
feature is missing in more metal-poor galaxies such as the 
LMC. Several theories have been put forward to explain the 
bimo dal distribution, but a satisfying scenario is still lacking 
(see lAntonello et aTl l2002. for an extended discussion). 

Since the MW and M 3 1 have similar morphological type 
and chemical composition, the observed period distributions 
are expected to be also similar. The results obtained so far 
have prevented a direct comparison because of observational 
biases, which are often difficult to evaluate. On the one hand, 
faint Cepheids (with usually short periods) can be missing in 
some shallow surveys and, on the other hand, the distribution 
of the observations can make the identification of long period 
Cepheids difficult. 

To e valuate the observational biases, the 416 Cepheids from 



PaperJ were compared with the 420 Cepheids in M31 of 



the GCVS dSamus et alj|2004 . Both samples have almost the 
same number of stars, but the two period distributions (Fig. 
[2]) are found to be largely different, as demonstrated by the 
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Fig. 2. Top: Norm alized pe riod distribution for the 416 de- 
tected Cepheids in IPaper B compared with the 420 Cepheids 
in M31 of the GCVS. Bottom: Normalized period distribu- 
tion for the 416 detected Cepheids in M 31 compared with the 
509 Cepheids in the MW from the David Dunlap Observatory 
(DDO) sample. The OGLE Large Magellanic Cloud Cepheid 
sample is also shown for comparison. 



Kolmogorov-Smirnov tests, which provide values of around 
10~ 13 . The origin of the observed difference can be explained 
by two main factors. Firstly, the Paper I sample is known to 
have a bias for the longest period Cephe ids because of the 
observational window function (see IPaper I for further infor- 
mation). And secondly, the GCVS presents an important bias 
at short periods because the GCVS is shallower than our sur- 
vey and faint Cepheids are missing. Furthermore, both Cepheid 
samples in M 3 1 can be compared with the MW samp le of the 



David Dunlap Observatory!] (DDO. lFernie etalJll995 



). As can 

be seen (Fig. [2j, the period distributions of IPaper II and the 



DDO sample are very similar (with a Kolmogorov-Smirnov 
test value of 0.49). Only a slight difference is observed at 
long periods, possibly because of the observational bias in the 
Paper J sample. Therefore, as expected, equivalent period dis- 
tributions are obtained for similar galaxies, including the sec- 
ondary peak at P > 10 d in metal rich galaxies. 

The conspicuous similarity between the period distribu- 
tions of MW and M 3 1 becomes even more striking when com- 
pared with a galaxy of different metallicity (Fig. [2]). More than 
1300 Cepheids ha ve been observed in LMC as part of the 
OGLE II surve>0 dUdalski et alJll999h . The resulting period 
distribution reveals a single peak which, as expected, is shifted 
towards shorter periods with respect to those of the distribu- 



1 Data obtained from: http : //www . astro . utoronto . ca/DD0/ 
research/cepheids/cepheids .html 

2 Data obtained from: ftp://sirius.astrouw.edu.pl/ogle/ 
ogle2/var_stars/lmc/cep/catalog/ 



tions of the MW an d M 3 1 . Therefore, the great similarity be- 
tween the DDO and lPaperll samples indicates that both the ob- 
servational biases and period distributions are equivalent. 

4. Fourier decomposition 

It is well known that Fourier decomposition of the Cepheid 
light curves can provide valuable information on the real na- 
ture of these variable stars, especially on their pulsation modes 
( Beaulieu et al.l[l995 ). The procedure involves fitting the coef- 
ficients (Ak and <pk) of a Fourier series of the form: 
./ 



m(t) - A + ^A k cos (27rKt>(f) + ip k ) 



(1) 



where m(t) and <t>(?) are the magnitude of the Cepheid and the 
phase at time t, respectively. Once the best coefficients have 
been obtained, the pulsation mode of the Cepheid variable can 
be derived from the amplitude ratio R^i - A k jA\ and phase 
difference ipt\ — fk — kipi- 

A prerequisite of the Fourier decomposition of the 416 
Cepheids identified in M31 is the determination of the 
maximum order of the fit, J. Simil a rly to other studies 
in the literature dAntonello et al.1 \l99(j iBeaulieu et al.l fl995 



Zakrzewski et al. 2000), we observed that an iterative solution 
provided satisfactory results. Therefore, for each light curve, 
the fitting procedure started with 7=1 and increased until 
the value of x\ t varied by less than one. Although the Fourier 
fit always provided values for all the light curves, several im- 
portant aspects were considered before studying the results. 
Firstly, since measurements with large photometric errors can 
provide unreliable fits, light curves with a mean photometric 
error larger than 0.1 mag were excluded. Secondly, all light 
curves W\i\\x do f > 7 were rejected to eliminate the fits that did 
not accurately match the observations. Finally, large gaps in a 
light curve can prevent the fit to correctly reproduce its shape, 
even when x\ t is relatively low. Problematic fits usually have 
coefficients with large error b ars. Therefore, coefficient errors 
were computed (according to Petersen 1986b and light curves 
with cr Ao > 0.01 mag, <x R21 > 0.1 or cr^,, > 0.6 rad (^ 0.27rrad) 
were also rejected. 

The above criteria yielded a sample of 315 Cepheids with 
accurate Fourier fits in both passbands (B and V). The R21 and 
(P21 values in both passbands (Pig. [3j have been used to classify 
75 Cepheids pulsating in the first-overtone mode (FO) and 240 
pulsating in the fundamental mode (FM). Since FO Cepheids 
are not expected to have periods longer than ~ 7 days, the FM 
sample has been completed with all Cepheids with log P > 0.9 
(^ 8 days), regardless of the quality of their Fourier fit. The fi- 
nal sample of 28 1 FM Cepheids includes three type II Cepheid 
candidates, but they are easily identified as a consequence of 
the detailed analysis of the P-L distribution (Sect. lBTTT i. 

5. Period-luminosity relationship 

The P-L diagram for the 356 classified Cepheids (28 1 FM and 
75 FO) reveals a large scatter in both passbands (Fig. 0J, es- 
pecially in B. The origin of the observed scatter can be un- 
derstood if we consider that the field of view covers about 7 
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Fig. 3. Fourier coefficients in each passband (B and V) as a 
function of period. The pulsation mode of the 315 Cepheids 
is identified: black circles - fundamental mode; grey circles - 
first overtone. 



kpc at the distan ce of M31 ((m - M) = 24.44 + 0.12 mag, 
Rib as et al 1l2005h . The large covered area in M 3 1 , the fact that 
Cepheids are located along the spiral arms, and the clumpy 
structure of the interstellar medium (clearly observed in the 
survey images), indicates that differential absorption in the disk 
is the most likely responsible for an important fraction of the 
observed scatter. In addition, the larger dispersion of the P-L 
distribution in B compared to V reinforces the interstellar ab- 
sorption hypothesis. 

Another possible effect on the observed scat ter is metal- 



Fig. 4. Observed magnitudes as a function of period for 356 
Cepheids. Filled black circles: Fundamental mode Cepheids 
with accurate Fourier coefficients. Empty black circles: 
Fundamental mode Cepheids without accurate Fourier coeffi- 
cients. Grey circles: First overtone mode pulsators. 



red giant branch and Cepheid distance determinations to sev- 
eral galaxi es has shown that a slight dependence on metallicity 
is present dSakai et al.l l2004). Therefore, a metallicity depen- 
dence on the zero point of the Cepheid P-L relationship could 
exist. Since metallicity withi n M 3 1 decreases as a function of 
the galactocentric distance (Zarit sky et al.lll994T) . some of the 
observed scatter could be as well introduced by metallicity dif- 
ferences. 

Also, the Cepheid measured photometry can be affected 
by blends. The mean seeing of our images is around 1 arc- 
sec. Therefore, each observed point-like source corresponds to 
~ 4 pc at the distance of M31. Since Cepheids are usually lo- 
cated in young star clusters and associations, it is likely that 
the observed magnitude of a Cepheid can be the combination 
of several unresolved sources. 

Below we analyse the three aforementioned effects that are 
the probable sources of the scatter in the P-L diagram. 

5.1. Absorption 

The effect of differential absorption can be partially corrected 
from the observed Cepheid color (B - V) through the color- 
excess (E(B - V) = (B - V) - (B - V) n . lFeasj [l999l). We 
used P-L relationships from Udalski et al. ( 1999H 3 ! to estimate 
Bq and Vo values. The LMC distance modulus was assumed to 
be (m - M)o = 18.42 + 0.06 mag and was computed from the 



licity. Several authors have suggested (see, i.e., iGieren et al 



2005) that the slope of the Cepheid P-L relationship is basically 
independent of metallicity. However, comparison of tip of the 



3 Updated relationships were obtained from the OGLE II web site: 
ftp : // sirius . astrouw . edu . pi/ ogle/ ogle2/var_stars/lmc/ 
cep/catalog/README . PL 
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Fig. 5. Absorption-corrected V magnitude as a function of pe- 
riod for 281 fundamental mode Cepheids. Filled circles: Stars 
with accurate Fourier fit. Empty circles: Stars without accu- 

1 19991) P-L relation- 



rate Fourier fits. Black line: Udalski et al 



ship, for a distance modulus to M 31 of 24.44 mag dRibas et al.l 
120051) . Grey lines: Dista nce corrected P-L r elationships for type 
II Cepheids as given bv lAlcock et al.l dl998l) . 



weighted mean of all the LMC distance determinations with 
eclipsing binaries (Fitzpatrick et al.ll2003l) . The reason for us- 
ing LMC relationships (instead of those for the MW), was 
motivated by several recent results. First. iGieren et al. (2005) 
suggested that the MW Cepheid distances, mostly obtained 
through the Baade-Wesselink method, could be affected by a 
systematic bias when converting th e observed radial veloci- 
ties into pulsation velocities. Second. lMacri et al.1 (120061) found 
that LMC P-L slopes provide a better agreement with the P- 
L distribution of Cepheids in NGC4258 than MW relation- 
ships. Fin ally, the new parallax measurements of severa l MW 
Cepheids (Be nedict et al . 2007; Ivan Leeuwen et al . 2007) yield 
P-L slopes that seem to be in better agreement with the LMC 
than with previous MW relationships. 

After obtaining the color excess for each Cepheid, a total- 
to-selective extinc tion ratio of %y = MV) /E(B - V) = 3 . 1 + 0.3 
dFitzpatrickll 19991) was used to compute the absorption and the 
Vo magnitude for every Cepheid. The resulting P-L distribution 
can be compared with the LMC P-L relationship (Fig. |5), as- 
sumi ng a M 3 1 distance modulus of (m - M)o = 24.44 + 0.12 
mag (Rib as et al. I l2005l) . An offset between the Vo values and 
the LMC P-L relationship is clearly observed, but the effects of 
metallicity and blending have still to be considered. 

Finally, as previously mentioned, three of the studied 
Cepheids seem to be of type II. To further investigate these ob- 
jects they have been kept in the studied sample of FM Cepheids 
(Sects. E3]and|6j. 



5.2. Metallicity 

As explained above, the zero point of the P-L relationship may 
depend on metallicity. To account for this effect, the metallicity 
of each Cepheid needs to be estimated, and we did so by con- 
sidering a galactoc entric metallicity dependence. According to 
Zari tsky et al. ( 1994 ). the M31 galactocentric metallicity de- 
pendence can be modeled by: 

12 + log(0/H) = (9.03 + 0.09) - (0.28 ± 0.10)(p/p - 0.4) 

where p is the de-projected galactocentric radius and po is the 
isophotal radius (77.44 arcm in for M31). Following the pro- 
cedure in iBaade & Arpl ( 1964b . the Cepheid galactocentric ra- 
dius was obtai ned assuming a po sition angle of 38°, an incli- 
nation of 1 2° 5 dSimien et al.ll 1978 ) and an M 3 1 center position 
of a = 00 /! 42 m 44 s .31 8 = +41°16'09'.'4 (in J2000.0 coordinates, 



Cotton et al.lll999h . 

The metallicity dependence of the Cepheid P-L relationship 
and the galactocentric metallicity variation implies that a differ- 
ent P-L relationship must be used for each individual Cepheid 
when computing the value Vq-M v . Alternatively, one can com- 
pute Vo - My by using a universal P-L relationship and then 
appl ying a correctio n for each Cepheid a posteriori. According 



to lSakai et al.l (120041) . the correction on the computed Vo - My 
value is 6(m-M)q/6(0/H) = -0.25 + 0.09 mag/dex when us- 
ing the lUdalski et al.l (1 19991) P-L relationships. Therefore, we 
computed metallicity corrections by comparing the assumed 
metallicity of each M31 Cephe id with the met allicity of the 
LMC (12 + log(0/H) = 8.5 dex. lSakai et alll2004 . The result- 
ing corrections for the studied sample lie between -0.15 mag 
and -0.05 mag. 



5.3. Blending 

The best approach to study the effect of blending given our 
available data set s is by using the observed Cepheid amplitudes 
( Antonellol l2002l) . The intrinsic amplitude M' of a Cepheid is 
given by the expression: 



M' = -2.5 log 



fc 



(2) 



where f c and F c are the fluxes at minimum and maximum light, 
respectively. When blending is present, the observed amplitude 
J\ can be expressed in the following form: 



m = -2.5 log 



fc+f 

Fc+f 



(3) 



where / is the sum of fluxes of all blending sources. From these 
two equations, one can write: 



J?l = 2.5 log 



1 + L 

1 + fc 



(4) 



Since M and f/f c are always positive quantities, the observed 
amplitude of a blended Cepheid is always smaller than its 
intrinsic amplitude. The observed amplitude and the blended 
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mean magnitude < to > can be expressed as: 



< m >= 



m - M 
m + M 



(5) 
(6) 



where to and M are the observed magnitudes at minimum and 
maximum light of a Cepheid, respectively. In the same way, 
the intrinsic amplitude and the blending-free mean magnitude 
< to' > can be expressed as: 



■■ m' - M l 
to' + M' 



< m >- 



(7) 
(8) 



From equations (HJl-® the difference on the mean magni- 
tude of a given Cepheid as a consequence of blending can be 
expressed as: 



A =< m > — < m' >— (to — to') — 



(9) 



Considering that (to - to') can be defined as: 

m - to'' = -2.5 log (^jA = ~ 2 -5 log |l + j J (10) 

and isolating f/f c from (0|, the difference on the mean magni- 
tude (0 can be expressed as: 



A = 2.5 log 



10' 



0.2SH 



- 10 



-0.2^1 



(ID 



Therefore, the variation on the mean magnitude of a given 
Cepheid because of blending can be computed from the ampli- 
tude of the Cepheid. Although M' and A are unknown quan- 
tities, the above equation can be solved in combination with 
period-color an d P-L relation s hips. The results shown below 
are based on the Udalski et al. (1999) relationships, but almost 
identical results wer e obtained with other LMC relationships 
dSandage et al.ll2004 . 

The next step of the process involves recalling that 
< m >- to, where to is the phase-weighted intensity-average 
mean magnitude of the observed Cepheids. Therefore, defin- 
ing the blending-free color excess as: 



£(fi' - V) s (B' - V) - (B - V) Q 

the blended color excess can be expressed as: 

E(B - V) = £(fi' - V) + A B - A v 



(12) 



(13) 



where A^ = B - B' and Ay = V - V. Analogously, the distance 
modulus can be expressed as: 



(Vo - M v ) 
where: 



(m - M)o + Ay - 7?y(A B - Ay) 



(Vo -M V ) = V- <R V E(B - V) 
(m - M) = T - KvEiB 1 - V) 



My 
My 



(14) 



(15) 
(16) 



The combination of Eqs. ( fT4b and ( fTTT i reveals that the value 
of (Vo - My) is only a function of the amplitude of the Cepheids 
and the distance modulus to M31. In addition, the fact that 
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Fig. 6. Amplitude relationship for the 240 fundamental mode 
Cepheids. 



the color excess is obtained from the observed color of the 
Cepheids (which is affected by blending) introduces the color 
term (A B - Ay) in Eq. ( TPfl i. Therefore, and contrary to the intu- 
itive interpretation, when reddening is computed from the ob- 
served color of Cepheids (B - V) the blended distance modulus 
can be either larger or smaller than the intrinsic distance mod- 
ulus (to - M)o, depending on the values of (A# - Ay). 

Considering that all Cepheids are r oughly at the sam e dis- 
tance, (to - M) = 24.44 + 0.12 mag dRibas et alj|2005b . and 
assuming that M' and J\' B are linearly dependent (Pig. [6j>, Eq. 
(fT4l can be numerically solved. In fact, when the intrinsic am- 
plitude of the Cepheid tends to zero, the observed amplitudes 
have to be zero in all passbands. Therefore, both amplitudes 
were considered to be proportional (i.e., ^{' B = aJ\' v ). The 
assumption that Cepheid amplitudes, when observed at dif- 
ferent passbands, are proportional can be considered a first- 
order approximatio n of the Fourier coefficient interrelations 
dNgeow et alJl2003h . 

Since the interrelations between B and V have not been ac- 
curately worked out, we computed the proportionality factor 
empirically from the observed amplitude of Cepheids. The am- 
plitude of the studied Cepheids can be reliably estimated from 
the Fourier fits (Sect.©. Therefore, the 240 FM Cepheids with 
accurate Fourier fits were used to compute the amplitude pro- 
portionality factor (Fig. |6j, obtaining a mean value (with 2.5<x 



clipping) of &1 B I9% WbI^v = 1-435 + 0.011. 

The above distance modulus to M31 and the amplitude 
proportionality factor, when applied to Eq. (fl4l . yielded the 
intrinsic amplitudes of the Cepheids. Once the intrinsic ampli- 
tudes are known, the difference on the mean magnitude can be 
computed (Fig. |7). The large resulting uncertainties in A# and 
Ay clearly indicate that not much can be said from the indi- 
vidual blending of each Cepheid. The uncertainties were ob- 
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Fig. 7. Observed amplitude of the 240 fundamental mode 
Cepheids with accurate Fourier fits as a function of the com- 
puted variation on the mean magnitude due to blending. Lines 
of constant intrinsic amplitude from 0.5 to 2 or 2.5 mag in steps 
of 0.5 mag are also shown (dashed lines). 

tained from a MonteCarlo run with 1000 realisations on the 
input parameters of each Cepheid. The positive A values (im- 
plying negative blending) provide additional information on the 
associated uncertainties (probably of the order of 0. 1-0.2 mag). 
In any case, it is interesting to observe that moderate intrinsic 
amplitudes are obtained for the entire sample, as it can be de- 
duced from the lines of constant intrinsic amplitude in Fig.|7](at 
A = 0, — M'). Furthermore, the Cepheid with Ay > 1 mag 
is the only suspected type II Cepheid with an accurate Fourier 
fit (see Fig. |5). The unrealistic blending value provides an ad- 
ditional evidence in favor of the type II classification. 

Given the large number of studied Cepheids, the obtained 
blending values can also be used to determine a mean blend- 
ing value for Cepheids in M31. The results are shown in 
Table [JJ where they are also compared with previous blend- 
ing determinations in M31 and M33. The values in iPaper I 
were obtained from t he third lig ht c ontribution in eclipsing bi- 
nary systems, and the MMSSOOl and lMMSS04l values were ob- 
tained from comparison of Hubble Space Telescope (HST) and 
ground ba sed images. T he blending factors (S) presented in 
MMSS00 and MMSS04 were transformed into variation on the 
mean magnitude values by considering that: 



A = -2.5 log 1 + 



/ 



<fc> 



= -2.5 log (1 +5) 



(17) 



where < f c > is the flux of the Cep heid at me an magnitude. 

The observed difference with MMSS00 could be due to 
the assumed distance moduli or to the observing conditions. 
On the one hand, the reported uncertainties on the LMC and 



M31 distance moduli could produce a variation of 0.12 mag 
in A B and Ay. On the other hand, blending depends on see- 
ing conditions and background level. Large seeing images in- 
crease the blending and high signal-to-noise data makes that 
faint stars are detected, decreasing the computed background 
and increasing the blending contribution. Considering that sim- 
ilar results have been obtained for the eclipsing binary sam- 
ple (which comes from the same observational data but from a 
completely different pr ocedure), the obtained differences with 
the values reported in MMSS00 could be the r esult of di ffer- 
ent observing conditions. The method used by|MMSS00 can 
only provide lower limits to blending values, since the HST 
point spread function (used as blending-free reference) could 
still hide unresolved companions. This is especially the case 
for the B-band results, where the available HST data was of low 
signal-to-noise, thus implying that only the most severe blends 
were dete cted. In addition, the relatively small sample used by 
MMSS00 cannot be considered representative of the Cepheid 



population in M 3 1 . 

It is interesting to note that the results of MMSS04, based 
on a relatively large sample of Cepheids in M 33, have mean 
values similar to ours (although slightly different median val- 
ues; TableHJ. Such good agreement should be taken with cau- 
tion as the host galaxies, the methods used, the period distribu- 
tions of the samples, and therefore the associated systematic er- 
rors are different. Nevertheless it is encouraging that estimates 
made using different methods in two somewhat differing spiral 
galaxies yield similar results. 

To further analyse the effects of blending we have applied 
several cuts to the FM sample. It has traditionally been ar- 
gued that the longer peri od (and brighter) C epheids should be 
less affected by blending (Macrietal . 2006). We computed the 
mean blending for Cepheids with a period longer than 12 days 
and obtained, as expected, lower blending values (Table [TJ. 
Furthermore, considering that blending decreases the Cepheid 
amplitude, larger amplitude Cepheids should also be less af- 
fected by blending and, in fact, the 37 Cepheids with J\y > 0.8 
mag do have lower blending values (Table [TJ. 

Finally, a blending-corrected color excess can be computed 
from Eq. (fTTt . obtaining a mean value of < E(B' - V) >= 
0.305 +0.01 1. The blended mean color excess is slightly lower 
(< E(B-V) >= 0.296+0.012), indicating that blending sources 
are bluer than Cepheids. The difference is more significant 
when considering only Cepheids with large blending values 
(Ay < -0.5 mag and A# < -0.5 mag), obtaining a blending- 
corrected color excess of < E{B' - V) >= 0.472 + 0.028, 
whereas < E(B - V) >= 0.315 + 0.030. Therefore, consid- 
ering the color-magnitude diagram (Fig.[TJ and the distribution 
of Cepheids along the spiral arms, it is likely that early-type 
and young main sequence stars are responsible for the large 
measured blendin g. The obta ined results are in good agreement 
with the results in MMSS04, who also found that M 33 blend- 
ing sources were on the average bluer than the Cepheids. To 
further check this scenario we used young MW open clusters 
to compute the light contribution that the main sequence stars 
would introduce on a Cepheid variable in the cluster. The re- 
sulting predictions on Ay and A# - Ay are in good agreement 
with the inferred blending values in M 3 1 . 
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Table 1. Blending values in M 31 and M 33 obtained from several methods. Top: Mean values and their errors. Bottom: Median 
values. 



Reference 


Sample 


Number 


< S v > 


< S B > 


< Ay > 


< A s > 












[mag] 


[mag] 


"\ MH x o o f\r\ 

MMSSOO 


Cepheids in M 3 1 


22 (10 in B) 


0.19 ± 0.03 


0.06 ± 0.06 


-0.18 ± 0.03 


-0.05 ± 0.05 


Paper I 


Eclipsing binaries 


48 


0.31 ± 0.07 


0.30 ± 0.06 


-0.23 ± 0.05 


-0.25 ± 0.04 


Present work 


FM Cepheids 


240 


0.31 ± 0.03 


0.37 ± 0.04 


-0.23 ± 0.02 


-0.24 ± 0.03 


Present work 


FM Cepheids with P > 12 days 


37 


0.14 ±0.06 


0.20 ±0.10 


-0.09 ± 0.05 


-0.10 ±0.07 


Present work 


FM Cepheids with !A V > 0.8 mag 


37 


0.10 ±0.05 


0.11 ±0.07 


-0.06 ± 0.04 


-0.05 ± 0.06 


MMS5U4 


t^epneicls in JV1 3d 


yj (j 1 in d) 


U./4 ± U.Ui 


U.2y ± U.Uo 


— U.ZU ± U.UZ 


— V.Zi ± U.U4 


MMSS04 


Cepheids in M 33 with P > 10 days 


60 (39 in B) 


0.16 ±0.04 


0.20 ± 0.05 


-0.14 ±0.025 


-0.17 ±0.04 


Reference 


Sample 


Number 


Median(5 v ) 


Median(5 fl ) 


Median(Av/) 


Median(A s ) 












[mag] 


[mag] 


MMSSOO 


Cepheids 


22 (10 in B) 


0.12 


0.00 


-0.12 


0.00 


Paper I 


Eclipsing binaries 


48 


0.09 


0.16 


-0.09 


-0.16 


Present work 


FM Cepheids 


240 


0.20 


0.15 


-0.19 


-0.16 


Present work 


FM Cepheids with P > 12 days 


37 


0.09 


0.04 


-0.09 


-0.05 


Present work 


FM Cepheids with &L V > 0.8 mag 


37 


0.09 


0.04 


-0.09 


-0.05 


MMSS04 


Cepheids in M 33 


95 (57 in B) 


0.13 


0.15 


-0.13 


-0.15 


MMSS04 


Cepheids in M 33 with P > 10 days 


60 (39 in B) 


0.07 


0.10 


-0.07 


-0.10 



6. Distance determination 

As previously mentioned, long period and large amplitude 
Cepheids are less affected by blends. Therefore, when com- 
puting the mean distance modulus to M31, a systematic trend 
should be observed for increasing period and amplitude cuts 
(i.e., rejecting short period or low amplitude Cepheids). The 
Vb _ My values obtained after the metallicity correction (Sect. 
15.2b were used to compute weighted mean distance determina- 
tions to M 3 1 (with 2.5cr clipping). The resulting values assume 
a distance modulus to the LMC of (m - M)o = 18.42 + 0.06 
mag (Sect. IBTTl i. The amplitudes for the 41 FM Cepheids with- 
out accurate fits were computed from the observations and also 
included in order to obtain a better coverage in periods and am- 
plitudes. 

Figure [8] shows that the distance modulus increases as the 
minimum period or minimum amplitude cuts increase. We ob- 
serve that the distance modulus value stabilizes for ^ly > 0.8 
mag. The most likely cause for the observed tendency is that 
blending is lower than photometric errors (or even zero) for 
large amplitudes. The most suitable period cut is more difficult 
to compute, although an increasing tendency is also observed. 
The observed behavior is explained if long period Cepheids are 
still affected by large blends, thus still introducing a bias in the 
derived distance modulus. Hence, an amplitude cut is the best 
choice to derive the distance to M 3 1 from the studied Cepheid 
sample. From this analysis (Fig. [8}, the 66 Cepheids with am- 
plitude J[ v > 0.8 mag seem to represent the best sample for 
distance determination. 

The effect of removing small amplitude Cepheids is clearly 
visible in Fig. [9] where the P-L diagram for the largest ampli- 
tude Cepheids is shown. Two stars with Jl v > 0.8 mag are 
placed far below the general P-L distribution (empty circles in 
Fig. |9). Both stars are on the type II Cepheids relationships, 
reinforcing the hypothesis that these stars are, in fact, type II 



24.4 



2^.35 



24.3 



24.25 



24.2 



24.1 



1.2 0.4 0.6 0.8 



1 1.2 



0.2 0.4 0.6 0.8 



Fig. 8. Weighted mean Vq - My values for the 28 1 fundamental 
mode Cepheids in M31 for different cuts. Error bars indicate 
the error on the mean. The dashed line represents the adopted 
distance determination of (m - M)o = 24.32 + 0.12 mag. Left: 
Each distance determination includes all Cepheids with period 
longer than the specified value. Right: Each distance determi- 
nation includes all Cepheids with amplitude larger than the 
specified value. 



Cepheids and they have not been considered from now on. 
Finally, the obtained P-L sl ope, with a va l ue of -2.83 + 0.12 
mag dex -1 , is closer to the lUdalski et alj (11999b LMC value 
of -2.779 + 0.031 mag dex" 1 than to MW P-L slopes (e.g.: 
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Fig. 9. Absorption-corrected V magnitude as a function of pe- 
riod for 281 fundamental mode Cepheids. Empty grey cir- 
cles: Complete fundamental mode sample. Filled black circles: 
Cepheids with 3Ky > 0.8 mag. Empty black circles: Two stars 
with ${ v > 0.8 mag excluded fr om the distance determination. 
Black line: ludalski et alJd 1999b P-L relationship, for a distance 
modulus to M 3 1 of 24.44 dRibas et al.l2005b and a mean metal- 
licity correction of -0.1 mag (Sect. lBT2"l l. Grey lines: P-L rela- 
tionships for type II Cepheids as given bv lAlcock et al. (1998) 
at distance of M 3 1 . 



-3.087 + 0.085 mag dex" 1 , Sanda ge et alj|2004b . therefore fa- 



voring the adopted P-L relationships. 

From the considerations above, the M 3 1 distance modulus 
obtained from the studied sample of Cepheids is (m - M)o — 
(Yo - Mv)m>Q.% = 24.32 + 0.12 mag. This value is compati- 
ble with most distance determinations found in the literature 
and the assumed distance modulus in Sect. 15.31 Considering 
that some large amplitude Cepheids may still be affected by 
blends, the distance modulus we derive could have a slight 
negative bias. Note that our final value is 0.09 mag larger than 
the weighted mean distance modulus for the 28 1 FM Cepheids 
((Vb _ My) - 24.23 +0.12 mag). Therefore, blending is clearly 
an important effect that has to be considered when obtaining ex- 
tragalactic distance determinations. In fact, blending has nearly 
the same impact on the final distance determination for M 3 1 as 
the metallicity correction. 

7. Summary 

As a resu lt of a va riability survey in the North-Eastern quadrant 
of M31 (Paper I), 416 Cepheids were detected and measured 
with a large number of epochs (~ 250) per filter (B and V). 
The large number of detected Cepheids has allowed the direct 
study of their period distribution. The results reveal that the 
period distributions in M 3 1 and the MW are extremely similar 



and, in addition, that the Cepheid sample obtained is almost as 
complete as the DDO MW sample. The only difference may be 
the lack of long period Cepheids. 

The large number of epochs obtained in both filters permits 
an accurate pulsation mode identification for 240 FM and 75 
FO Cepheids. Althoug h some FO Cepheids were previously 
detected in M31 (IFliri et alJl2006h . our sample represents an 
important increase in the number of detected FO pulsators and 
opens a new window to study the basic properties of these stars 
in another Local Group galaxy. The sample of FM Cepheids 
was completed with 41 long period (i.e., longer than 8 days) 
Cepheid variables, resulting in 28 1 FM Cepheids. The posterior 
analysis revealed that at least three of these stars are type II 
Cepheids. 

The analysis of the P-L relationship for the FM Cepheids 
reveals a large scatter, which is not explained solely through 
the effects of interstellar absorption and metallicity. Although 
additional efforts are needed to reduce the obtained uncertain- 
ties, a new method to compute the effect of blending is pre- 
sented. The exact dependence of the amplitudes among differ- 
ent passbands, as well as more precise amplitude determina- 
tions, would greatly improve the results shown here. Even by 
considering the large uncertainties, the large number of stud- 
ied Cepheids provides an accurate characterization of the mean 
blending values. We conclude that the most likely cause of 
blending seems to be the light from unresolved stars belonging 
to the same stellar associations or clusters as Cepheid variables. 

The effect of blending has been shown to be larger than 0.09 
mag in the distance modulus to M31, thus having an effect 
as important as the metallicity correction. Therefore, blend- 
ing should always be taken into account when obtaining ex- 
tragalactic distance determinations with Cepheids. It has also 
been shown that an amplitude cut is more adequate to obtain- 
ing unbiased distance determinations than a period cut, since 
the effect of blending is better accounted for. Finally, the large 
amplitude Cepheids have been used to obtain a distance modu- 
lus to M31 of (m - M) = 24.32 + 0.12 mag. 
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